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Monopole and photon contributions to Wilson loops are calculated using Monte-Carlo simula- 
tions of SU(2) QCD in the maximally abelian gauge. The string tensions of SU(2) QCD are well 
reproduced by extended monopole contributions alone. 
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To understand color confinement mechanism in QCD 
is absolutely necessary for us to analytically explain 
hadron physics starting from QCD The 'tHooft 

idea of abelian projection in which a partial gauge-fixing 
Q«.^ I is done keeping the maximal abelian torus group unbro- 
I ' ken is very interesting ||]. Then QCD can be regarded 
^ ^ asaC/(l)x[/(l) abelian gauge theory with magnetic 
O ■ monopoles and electric charges. An interesting gauge 
has been found among infinite ways of gauge-fixing for 
the abelian projection. It is called maximally abelian 
(MA) gauge in which link gauge fields are forced to 
become abelian as much as possible. 
PsJ ' In the preceding note, the authors js) have shown in the 
K*" i MA gauge and in SU{2) QCD that entropy dominance 
■ over energy of the monopole loops, i.e., condensation of 
the monopole loops occurs in the confinement phase if ex- 
tended monopoles are considered. After the abelian 
projection in the MA gauge, infrared behaviors of SU (2) 
QCD seem to be described by a compact-QED like U{1) 
theory with the running coupling constant instead of the 
, bare one and with the monopole mass on a dual lattice. 
^ • The confinement in SU{2) QCD may be interpreted as 
I ' the (dual) Meissner effect due to the abelian monopole 
condensation. 

If the monopoles alone are responsible for the confine- 
ment mechanism, the string tension which is a key quan- 
tity of confinement must be explained by monopole con- 
tributions. This is realized in compact QED as shown 
^ ■ recently by Stack and Wenslcy . The aim of this note 
5^ \ is to show that the same thing happens also in SU{2) 
QCD by means of evaluating monopole and photon con- 
tributions to Wilson loops. 

After the abelian projection in the MA gauge, a diag- 
onal matrix u{s,n) can be extracted uniquely from the 
original SU{2) link field. The diagonal matrix m(s,/j,) 
corresponds to a U{1) gauge field written by an angle 
variable 0^(s). 

Now we show an abelian Wilson loop operator (which 
we consider after the abelian projection) is rewritten by 
a product of monopole and photon contributions. Here 
we take into account only a simple Wilson loop, say, of 
size I X J. Then such an abelian Wilson loop operator is 
expressed as 

PF = exp{z^J,,(s)0^(s)}, (1) 
where J^{s) is an external current taking ±1 along the 
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Wilson loop. Since J^(s) is conserved, it is rewritten for 
such a simple Wilson loop in terms of an antisymmet- 
ric variable M^i^(s) as J„{s) = 9^',Af^^(s), where d' is a 
backward derivative on a lattice. M^,y(s) takes ±1 on a 
surface with the Wilson loop boundary. Although we can 
choose any surface of such a type, we adopt a minimal 
surface here. We get 

W^exp{-'-Y,M,As)UAs)}, (2) 

where ffi,y{s) = d^6^{s) — 8,^0 ^{s) and is a forward 
derivative on a lattice. The gauge plaquette variable 
can be decomposed into f^iy{s) = f^iy{s) + 2TTn^^(s) 
where f^v{s) G [— tt, tt] corresponds to a field strength 
and n^v{s) is an integer- valued plaquette variable |H| de- 
noting the Dirac string. Since M^,j{s) and na,j{s) are 
integers, the latter does not contribute to Eq. (0). Hence 
fij.v{s) in Eq. (^) is replaced by f^u{s). Using a decom- 
position rule 

we get 

W ^Wi-W2 (3) 
= exp{-i^9;/^,(s)i?(s - s')J.(s')} 

where a monopole current fc^(s) is defined as k^{s) — 
(l/47r)e^Q^^9ct/^-y(s) following DeGrand-Toussaint [[l2[ . 
D{s) is the lattice Coulomb propagator. Since f^i.u{s) cor- 
responds to the field strength of the photon field, VFi(W^2) 
is the photon (the monopole) contribution to the Wilson 
loop. 

To study the features of both contributions, we evalu- 
ate the expectation values (Wi) and (W2) separately and 
compare them with those of W . 

The Monte-Carlo simulations were done on 24!^ lattice 
from (3 = 2.4 to (3 = 2.8. All measurements were done ev- 
ery 30 sweeps after a thermalization of 1500 sweeps. We 
took 50 configurations totally for measurements. The 
gauge-fixing criterion is the same as done in Ref. [ p^ . 
Using gauge-fixed configurations, we evaluated monopole 
currents. As shown in the previous note type-2 ex- 
tended monopole loops with 6 > 6c ~ 5.2 x IO^'^(Al)^^ 
condense, where b = na{(3) for n"^ extended monopoles 
and a{(3) is the lattice constant. So we measured 2^ ex- 
tended monopole with 6 = 2a{(3) of the type-2 [0. Then 
the effective (renormalized) lattice volume becomes 12^. 
Since the original lattice is 24"*, 2'^ extended monopoles 
are the largest from which we can get useful data of the 
static potentials from Wilson loops. For (3 — 2.7 and 
2.8, the value b = 2a{f3) becomes less than be and so the 
monopoles may not reproduce the string tension. 
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We have evaluated the averages of W using abehan hnk 
variables (called abehan), of Wi ■ W2 (called total), of Wi 
(photon part), and W2 (monopole part), separately. Both 
the first and the second averages are evaluated to check 
reliability of the data, since both should be equivalent as 
known from Eq. (^). 

The results are summarized as follows. 

1. The monopole contributions to Wilson loops are 
obtained with relatively small errors. Surprisingly 
enough, the Creutz ratios of the monopole contri- 
butions are almost independent of the loop size as 
shown partially in Table |. This means that the 
monopole contributions are composed only of an 
area, a perimeter, and a constant terms. 

2. Assuming the static potential is given by a linear 
+ Coulomb + constant terms, we can determine 
them from the least square fit to the Wilson loops 
0. We plot their data in Fig. |l| (at /3 = 2.5) 
and in Fig. ^ (at /3 = 2.6). We find the monopole 
contributions are responsible for the linear-rising 
behaviors. The photon part contributes only to the 
short-ranged region. There seems to exist a small 
discrepancy between the abelian and the monopole 
-I- photon parts for R/a = 12, but finite-size effects 
are expected there. Similar data are obtained for 
/? = 2.4. 

3. This is seen more clearly from the data of the string 
tensions which are determined from the static po- 
tentials. They are shown in Fig. ^. Systematic 
errors coming from various least square fits are not 
completely certain and are not plotted in the figure, 
although they are not negligible. The string ten- 
sions are well reproduced by the monopoles alone 
for /3 < 2.6 and the photon part has almost vanish- 
ing string tensions. The string tensions evaluated 
from the total part (which are not shown here) are 
consistent with those of abelian and monopoles. 

A. At f3 — 2.7 and 2.8, the monopoles which have 
6 < 6c do not seem to reproduce the abelian string 
tensions as shown in Table However the string 
tensions from the total part which should be equal 
to the abelian ones are also smaller. The origin 
of the difference resides in the smallness of the 
renormalized lattice volume. To check it, we have 
adopted only even-sized abelian Wilson loops which 
correspond to the total ones with the lattice spac- 
ing 2a(/3) and made the least-square fit. The string 
tensions are almost unchanged for /? = 2.4, 2.5, 
and 2.6 (see Fig. ||), but they become smaller for 
P = 2.7 and 2.8 with large errors. The difference 
between the abelian (even only) and the total be- 
comes smaller. In conclusion, to get definite results 
for /3 > 2.7 in this framework of the analysis, we 
have to adopt much larger lattices like 48^* . We can 
not conclude at present that the monopoles with 
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b > be alone can reproduce the string tensions as is 
expected. 



5. We have derived also Coulomb coefficients from 
the static potentials as shown in Fig. ^. The 
monopole part has almost vanishing Coulomb co- 
efficients which is in agreement with the constant 
behaviors of the Creutz rations of the monopole 
part as shown above. The photon part has large 
coefficients, but they do not reproduce the coef- 
ficients of the abelian static potentials. But again 
they are not far from the Coulomb coefhcients from 
Eq. (H) (called total in Fig. ||) and those from even- 
sized abelian Wilson loops. The following may be 
interesting. The photon parts are evaluated on an 
effective lattice with b = 2a{(3). Hence they have 
different values of b for differnt (3. The Coulomb 
coefficients of the photon parts are well reproduced 
by the SU{2) running coupling constants g{b) with 
b = 2a(/?), i.e., -g{b)'^ /IGn, where 



The scale parameter A determined is A 46Al 
which is quite near the value A ~ 42Al fixed from 
the monopole action. 

Our analyses in this note as well as in the previ- 
ous one show that abelian monopoles are responsible 
for confinement in SU (2) QCD and condensation of the 
monopoles is the confinement mechanism if the abelian 
projection is done in the MA gauge. 

Finally we make some comments on the results of both 
notes. 

1. Why is the MA gauge so nice? Note that an abelian 
projection reduces QCD into an abelian theory 
with diagonal gluons as a photon-like particle and 
off-diagonal gluons as charged partices. The char- 
acteristic features of the MA gauge among many 
abelian projections are in the following. The MA 
gauge is defined in such a way that a quantity 



is minimized, where Ui{s,fi) and 6^2(5, m) com- 
ponents of a SU{2) hnk field U{s,fi) = Uo{s,fi) + 
U{s,fi) ■ a. Hence the gauge condition forces as 
many link fields as possible to become abelian, i.e., 
diagonal. However one can not let all link fields 
diagonal, i.e., R = as seen from a histogram anal- 
ysis shown in Ref. jl^ and from the correlation of 
R with the monopole density [^sj. If we see only 
long-ranged physics, SU{2) QCD in the MA gauge 
may be regarded as a [/(I) theory with abelian link 
photon fields and abelian monopoles. 
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Is there any other gauge showing similar behaviors? 
There seem to exist many candidates. For example, 
the following gauge may be interesting in which 

is minimized. The condition also forces link 
fields to be diagonal as much as possible. The 
case n — 2 leads us to J2s fj,^i(^)^Ji(''')i^t^ ^ 
igA^{s))A^{s) = in the continuum limit. The 
study in the gauge is in progress. 

2. Using the method developed in Refs. C^Jl^, the 
abelian monopole action derived in may be 
mapped on to a field theoretical model of an abelian 
(dual) Higgs system. The model is just equal to a 
Ginzburg-Landau type theory which one of the au- 
thors (T.S.) derived earlier assuming abelian dom- 
inance and monopole condensation p^ . 

3. To extend our method to a T 7^ system and also 
to SU{3) with or without dynamical quarks is very 
interesting. These studies are also in progress. 

We wish to acknowledge Yoshimi Matsubara for use- 
ful discussions especially on the entropy decreasing ef- 
fects of the monopole loops. This work is financially 
supported by JSPS Grant-in Aid for Scientific Research 
(c)(No.04640289). 



* e-mail addressishibaOhep.s. kanazawa-u.ac.jp 
** e-mail address:suzuki@hep.s. kanazawa-u.ac.jp 
[1] G.'tHooft, in High Energy Physics, edited by A.Zichichi 

(Editorice Compositori, Bologna, 1975). 
[2] S. Mandelstam, Phys. Rep. 23C, 245 (1976). 
[3] G. 'tHooft, Nucl. Phys. B190, 455 (1981). 
[4] A.S. Kronfeld et al., Phys. Lett. 198B, 516 (1987); A.S. 

Kronfeld et al., Nucl.Phys. B293, 461 (1987). 
[5] T. Suzuki and I. Yotsuyanagi, Phys. Rev. D42, 4257 

(1990). 

[6] S. Hioki et al., Phys. Lett. 285B, 343 (1992). 

[7] T. Suzuki, Nucl. Phys. B(Proc. Suppl.) 30, 176 (1993) 

and see also references therein. 
[8] H.Shiba and T.Suzuki, Kanazawa University, Report No. 

Kanazawa 93-09, 1993. 
[9] This condition f^v{s) £ [— 7r,7r] is applicable only in the 

case of defining a smallest 1^ monopole. 
[10] T.L. Ivanenko et al., Phys. Lett. 252B, 631 (1990). 
[11] J.D. Stack and R.J. Wensley, Nucl. Phys. B371, 597 

(1992). 

[12] T.A. DeGrand and D. Toussaint, Phys. Rev. D22, 2478 
(1980). 

[13] S. Hioki et al.,Phys.Lett. 271B, 201 (1991). 
[14] S.Itoh et al., Phys. Rev. D33, 1806 (1986). 



5 



[15] T.Suzuki and I.Yotsuyanagi, Nucl. Phys. B(Proc. Suppl.) 

20, 236 (1991). 
[16] S. Samuel, Nucl. Phys. B154, 62 (1979). 
[17] M. Stone and P.R. Thomas, Phys. Rev. Lett. 41, 351 

(1978). 

[18] T. Suzuki, Prog. Thcor. Phys. 81, 929 (1988); S. Maedan 
and T. Suzuki, Prog. Theor. Phys. 81, 229 (1989); T. 
Suzuki, Prog. Theor. Phys. 81, 752 (1989); S. Maedan 
et al, Prog. Thcor. Phys. 84, 130 (1990); S. Kamizawa 
et al., Nucl Phys. B389, 563 (1993); H. Monden et al., 
Phys. Lett. 294B, 100 (1992). 

[19] J. Smit and A.J. van der Sijs, Nucl. Phys. B355, 603 
(1991); Nucl. Phys. B(Proc. Suppl.) 20, 221 (1991). 



TABLE I. Croutz ratios from abchan and monopolc Wilson 
loops at /3 = 2.6. The Monopole Creutz ratio values are 
devided by 4, being adjusted to those in unit a(/3). 



Creutz ratios 


abelian 


monopole 


X(2,2) 


0.0872(2) 




X(3,3) 


0.0460(3) 




X(4,4) 


0.0323(4) 


0.0173(2) 


X(5,5) 


0.0270(7) 




X(6,6) 


0.0248(9) 


0.0169(3) 


X(7,7) 


0.0255(13) 




X(8,8) 


0.0196(25) 


0.0172(4) 


X(9,9) 


0.0147(66) 




X(10, 10) 


0.0307(160) 


0.0178(6) 


x(ii,ii) 


undeterminable 




X(12, 12) 


undeterminable 


0.0178(9) 



TABLE IL String tensions cr/Ai at /3 = 2.7 and 2.8. 
Abelian and total mean those from Wilson loops evaluated us- 
ing usual link fields and Eq. (3) respectively. Abelian (even) 
means the least-square fit using only even-sized abelian Wil- 
son loops. 



/? 


abelian 


abelian (even) 


total 


monopole 


photon 


2.7 
2.8 


1962(37) 
2130(32) 


1467(220) 
1621(187) 


1238(516) 
1269(326) 


1087(394) 
1025(261) 


-36(62) 
-21(134) 



FIG. 1. Static potentials aV{R) versus R/a at f3 = 2.5. 
The values are shifted by a constant. 

FIG. 2. Static potentials aV{R) versus R/a at /3 = 2.6. 
The values are shifted by a constant. 



FIG. 3. String tensions at /3 = 2.4, 2.5, and 2.6. 



FIG. 4. Coulomb coefficients. 
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